Quasielastic neutron scattering of hydrated BaZr .9ov4o.io02.95 (A = Y and Sc) 
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Proton motions in hydrated proton conducting perovskites BaZro.9v4o.io02. 95 (A = Y and Sc) have 
been investigated using quasielastic neutron scattering. The results reveal a localized motion on the 
ps time scale and with an activation energy of ~10-30 meV, in both materials. The temperature 
dependence of the total mean square displacement of the protons suggests an onset of this motion 
at a temperature of about 300 K. Comparison of the QENS results to density functional theory 
calculations suggests that for both materials this motion can be ascribed to intra-octahedral proton 
transfers occurring close to a dopant atom. The low activation energy, more than ten times lower 
than the activation energy for the macroscopic proton conductivity, suggests that this motion is not 
the rate-limiting process for the long-range proton diffusion, i.e. it is not linked to the two materials 
significantly different proton conductivities. 



I. INTRODUCTION 

Hydrated acceptor-doped perovskites are known to 
be proton conductors in the temperature range ~200- 
700° CP The doping creates an oxygen-deficient struc- 
ture and in a humid atmosphere water molecules dissoci- 
ate into hydroxide ions, which fill the oxygen vacancies, 
and protons, which bond to lattice oxygens P The pro- 
ton conductivity of perovskites varies quite dramatically 
amongst different types of materials for the same dopant 
concentration. For instance, in the BaZro.go^4o.io02.95 
(A = Y, In, Sc and Ga) system, the proton conductiv- 
ity differs several orders of magnitude depending on the 
dopant atom usedP The origin for this difference is in- 
triguing and remains to be clarified. 

The current understanding of the conduction mecha- 
nism in hydrated perovskites is that the transport pro- 
cess occurs in two elementary steps: i) hydrogen-bond 
mediated proton transfer between adjacent oxygens, and 
ii) reorientational motion of the hydroxyl group in be- 
tween such transfers. In this two-stage mechanism it is 
only the protons that exhibit long-range diffusion while 
the oxygens reside on their crystallographic positions. 
These processes have been studied in differe nt syst ems 
mainly using molecular dynamics simu lati on d 3 * 4 * 5 ^ and 
quasielastic neutron scattering (QENS^P^. QENS is a 
particularly suitable experimental technique as it gives 
access to the time range ^10 -13 -10 -9 s, in which these 
processes occur, as well as providing information about 
the spatial geometry of the atomic motions. In addition, 
the very high neutron scattering cross section of protons 
provides a good contrast in the experiments and enables 
studies of systems with a low dopant concentration, i.e. 
a low hydration degree and concomitantly few protons. 



Previous QENS investigations have revealed localized 
motions in the ps time range, mainly interpreted as ro- 
tational motions with radii between r~ 0.7 and 1.1 A, 
and with activation energies E a <100 meVP^The results 
have also been interpreted in more complex terms of a dif- 
fusional process of the proton consisting of a sequence of 
free diffusion and trapping/escape events PI These stud- 
ies have, however, been performed on orthorhombic per- 
ovskites, which makes the interpretation of the data more 
difficult as the structure is intrinsically anisotropic and 
there are several non-equivalent oxygen sites for the pro- 
tons in the structure. In addition, no study has addressed 
the question of the role of the type of dopant on the ele- 
mentary steps in the conduction mechanism. 

Molecular dynamics simulations on proton conducting 
perovskites have revealed a fast rotational motion of the 
protons around the line connecting the cell edge of the 
perovskite structure on the ps tim e scale and with an 
activation energy less than 50 meVPEEm By extending 
the simulations to higher temperatures and pressures, 
also five proton transfer events within a simulation time 
of ^100 ps were observedP This process was found to 
be considerably dependent on the oxygen-oxygen separa- 
tion, inferred by the vibrations of the oxygen sublatticeP 
In the contracted transition state, having an energy of 
~0.41 eV, proton transfer was found to occur almost bar- 
rierless, although it was not always that the proton was 
transferred.^ 

In the present work we investigate the proton dynam- 
ics in 10% Y- and Sc-doped BaZrC>3, the two extremes 
in the cubic BaZro.gov4o.io02. 95 (A = Y, In, Sc and Ga) 
system considering the conductivity, using QENS. The 
main interest in this study lies in the fact that although 
the structure of these materials is very similar, their pro- 
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ton conductivities differ almost two orders of magnitude 
with the Y-doped perovskite exhibiting the highest!^ The 
QENS results reveal a localized proton motion, with a 
slightly different activation energy in the two materials, 
and a comparison to density functional theory calcula- 
tions suggests that for both materials we observe proton 
transfers between neighboring oxygens close to dopant 
atoms. 



II. EXPERIMENTAL DETAILS 



The BaZr .9o^4o.io02.95 (^4 = Y and Sc) samples, here- 
after abbreviated as 10Y:BZO and fOSc:BZO, respec- 
tively, were prepared by mixing stoichiometric amounts 
of BaCC>3, ZrC>2 and (Y/Sc)203. The oxides were heated 
to 250°C overnight to remove moisture prior to weighing. 
To ensure thorough mixing ethanol (99.5 %) was added 
during the milling procedure, which was performed man- 
ually using an agate mortar and a pestle. The finely 
ground mixtures were fired at f000°C for 8 hours and 
subsequently ground and pelletized using a f3 mm di- 
ameter die under a pressure of 8 tons. The pellets were 
sintered at 1200°C in air for 72 hours. After sintering, the 
pellets were reground, compacted, and retired at 1500°C 
for 48 h. Finally, the pellets were finely reground to pow- 
ders. The charging with protons was performed by an- 
nealing the powder samples at ~300°C under a flow (12 
ml/min) of Ar saturated with water vapor at 76° C for 10 
days. Our X-ray measurements revealed a cubic struc- 
ture of both hydrated materials. One should note that 
10Y:BZO has previously been reported as tetragonalP 
However, from our diffractogram we find no evidence for 
a tetragonal distortion. 

The QENS experiments were performed at the time-of 
flight spectrometer IN6 at the Institut Laue Langevin in 
Grenoble (France), using neutrons with a wavelength of 
5.1 A, yielding an energy resolution of 100 /zeV (FWHM) 
and a Q-range 0.26-1.9 A- 1 . The Q-dependent spectra 
were summed into 14 groups to increase the statistics. 
Spectra were measured at 2, 250, 300, 380, 430, 465 and 
495 K for 10Y:BZO, and at 2, 380, 430, 465 and 495 K 
for 10Sc:BZO, and the integration time was four hours. 
The 2 K measurements were used as resolution functions 
in the data analysis. The samples were loaded in vacuum 
tight Al-containers, which were coated on the inside with 
a 100 nm thick layer of Pt to avoid corrosion. The sam- 
ple thickness was chosen to 6 mm. The sample containers 
were mounted in 135° geometry in a cryofurnace. A spec- 
trum of the empty sample container was measured at 380 
K and subtracted as a background from the sample spec- 
tra while a vanadium standard was used to correct for 
the detector efficiency. 



III. THEORETICAL BACKGROUND 

In a QENS experiment one obtains the dynamic struc- 
ture factor, S(Q,uj), which gives the probability that an 
incident neutron is scattered by the sample with a mo- 
mentum transfer HQ and an energy transfer fko. In our 
analysis we have modeled S(Q, u>) with the following scat- 
tering function 



Sq(Q,tj) = b(Q)-5(hw) 
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Here, the first term describes the elastic scattering from 
those atoms that move too slow to be resolved in the 
experiment {b{Q) is a Q-dependent constant and 5{hu) 
is a delta function). The second term, which appears as 
a broadening of the elastic peak, accounts for stochastic 
motions (translational and rcorientational) . This result 
in a small energy loss/gain of the incident neutrons. At 
the IN6 spectrometer motions in the time range ^10 -13 
to 10 -10 s can be studied. Since only protons are ex- 
pected to move on this time scale the quasielastic broad- 
ening can thus be directly related to proton motions, 
which facilitates the analysis. The quasielastic scatter- 
ing is described by a Lorentzian with amplitude a\ J (Q) 
and width Alu(Q) (HWHM) while c(Q) and d(Q) are pa- 
rameters of a sloping background. It should be noted 
that all fit parameters were free throughout the fitting 
procedure. 

The measured scattering function, S meas (Q, uj), is a 
convolution of the real scattering function, S(Q,uj), and 
the resolution function of the instrument, R(Q,lo), i.e. 



S mcas (<2,^) = S(Q,u)®R(Q,w). 



(2) 



Thus, the modeled scattering function in Eq. ([I]) is con- 
voluted with the resolution function in the data analysis. 

The Q-dependences of Alu and <2l contain informa- 
tion about the relaxation time and spatial geometry of 
the proton dynamics. For instance, for a long-range dif- 
fusional process and for sufficiently small Q- values, the 
quasielastic width follows a characteristic Q 2 -dependcncc 
according to Auj = hDQ 2 , where D is the diffusion 
constant.^ For a local process, e.g. a rotational mo- 
tion, Aw is practically Q-independent and is related to 
the relaxation time r through r = hj Aw One should 
here note that before analyzing the Q-dependence of the 
quasielastic intensities, ol, these need to be corrected 
for the overall decrease in scattering intensity due to the 
Debye- Waller factor, that accounts for the harmonic (vi- 
brational) motions. The total scattering can be written 
as 



S(Q,uj) = e~ 



~S q (Q,u>) 



(3) 



where (it 2 )h arm is the harmonic mean square displace- 
ment of all atoms in the material. Thus, also the mean 



3 



Quasielastic part 



Instrumental 
resolution 




-0.5 0.5 

co/meV 



FIG. 1: S meaa (Q,u) of 10Y:BZO and 10Sc:BZO at Q=1.9 
A" 1 and for T = 380-495 K. The solid lines represent fits 
to the experimental data (markers) according to Eq. Q. 
Each spectrum has been normalized to unity and cut at 
S m eas(Q,w) = 0.15 for easier comparison. 



square displacement is needed for a full analysis of the 
experimental data. Following a custom route in the anal- 
ysis of QENS data, the total mean square displacement 
(u 2 ), can be calculated from the ratio of the elastic and 
total scattering intensity according to 
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IV. RESULTS 



)Q 2 



(4) 



Figure [T] shows the experimental data, S mcas (Q , uj) , for 
the two investigated materials for temperatures in the 
range 380-495 K, and for Q = 1.9 A -1 . As seen in the 
figure there is a quasielastic component together with 
a strong elastic component for both materials in each 
spectrum. The quasielastic component increases with in- 
creasing temperature. 

To determine when we can actually expect to observe 
the stochastic motion of the protons in the experiment 
we have calculated the total mean square displacement, 
(it 2 ), using Eq. Figure [2] shows the logarithm of 
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FIG. 2: Q 2 -dependence of the logarithm of the elastic inten- 
sity normalized to the total intensity, where the error bars of 
the data points are within the size of the symbols. The lines 
represent linear fits with a slope -{u 2 }. 



the ratio of the elastic and total intensity versus Q 2 , 
together with linear fits from which (u 2 ) has been de- 
termined. We have discarded the spectra for 1.35 A -1 
< Q < 1.75 A -1 due to a Bragg reflection resulting in 
a strong elastic intensity masking the quasielastic con- 
tribution. The total intensity is approximated by inte- 
grating the spectrum over the region -2 < lu < 2 meV. 
Figure [3] shows the temperature dependence of the so ob- 
tained (u 2 ). For 10Sc:BZO we have only performed ex- 
periments at 2 K and T > 380 K, whereas for 10Y:BZO, 
we have also data for two temperatures in between. For 
10Y:BZO at low temperatures, T < 300 K, the mean 
square displacement increases only slightly and basically 
linearly with temperature, while above 300 K the temper- 
ature dependence is much stronger. The linear increase 
at low temperature reflects an increasing amplitude of 
harmonic vibrations while the non-linear increase above 
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FIG. 3: Temperature dependence of the total mean square 
displacement {u 2 }. The dashed lines represent the estimated 
harmonic mean square displacement, (w 2 )harm, for the two 
materials. 



300 K results from that we also resolve proton motions 
other than vibrational. The same behaviour can be in- 
ferred for 10Sc:BZO, even though we here do not observe 
the transition in the mean square displacement from a 
vibrational to a stochastic motion due to the lack of data 
in this region. Over the whole temperature range, the 
mean square displacement is slightly larger in the Sc- 
doped material compared to the Y-doped counterpart. 

Consistent with the mean square displacement we in- 
deed observe a quasielastic broadening for T > 300 K. 
Figure [4] shows the quasielastic widths as a function of 
Q, obtained from the fits of Eq. ([T| to the spectra. As 
seen, the widths are practically Q-independent and have 
values in the range 0.17-0.21 meV. There is a slight tem- 
perature dependence of the quasielastic widths. In or- 
der to extract this dependence we have fitted the largely 
Q-indcpcndent widths by a constant at each tempera- 
ture, represented by the lines in Fig. [4] An Arrhenius 
plot of the so obtained quasielastic widths are shown in 
Fig. [5] From this figure we obtain an activation energy 
of approximately 10 meV for 10Y:BZO and 30 meV for 
10Sc:BZO. As also seen in the figure, the width is larger 
for the Sc-doped material for all temperatures except for 
T = 380 K where they are the same. 

Figure [6] shows the Debye- Waller corrected quasielas- 
tic intensities versus Q for T = 380-495 K, for the two 
materials. It should be noted that the total mean square 
displacement in Fig. [3] cannot be used directly to calcu- 
late the Debye- Waller factor for these temperatures, as 
this should only contain the harmonic, vibrational, con- 
tribution. Therefore, for T — 380-495 K we have calcu- 
lated the Debye- Waller factor from a linear extrapolation 
of the vibrational mean square displacement between 2 
and 250 K for the 10Y:BZO material. For 10Sc:BZO. we 
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FIG. 4: Plots of the quasielastic widths and corresponding 
linear fits at 380 K (o), 430 K (A), 465 K (0) and 495 K 
(□). The data points and fits have been vertically shifted by 
0.1 meV (430 K), 0.2 meV (465 K) and 0.3 meV (495 K) for 
clarity. 



have estimated the harmonic mean square displacement 
by using the same slope as found for 10Y:BZO and an 
extrapolation from the 2K value, see Fig. [3] As seen in 
Fig. [6] the quasielastic intensity increases with increasing 
temperature and with Q. We find that the quasielastic 
intensity is slightly higher for the Sc-doped material. 



V. DISCUSSION 

The Q-independence of the quasielastic width observed 
in our experiment for both 10Y:BZO and 10Sc:BZO sug- 
gests that we arc following a local proton motion. 10 Based 
on the values of the widths (0.17-0.21 meV) this motion 
has a relaxation time in the range 3.1-3.9 ps, depend- 
ing on the temperature. The temperature dependence 
of the total mean square displacement in Fig. [3] suggests 
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FIG. 5: Arrhenius plot of the quasielastic widths with lin- 
ear fits which correspond to activation energies of 10 meV 
(10Y:BZO) and 30 meV (10Sc:BZO), respectively. 



an onset of this motion at a temperature of about 300 
K, whereas below this "threshold" temperature, the to- 
tal mean square displacement seems to be related only 
to harmonic (vibrational) motions of the protons. This 
finding agrees well with a recent inelastic neutron scatter- 
ing study of hydrated BaIno.20Zro.soO2.90, where it was 
found that the total mean square displacement of the 
protons increased only very slightly between 30 and 300 

kP 

The observed local process has an activation energy 
of 10 meV (10Y:BZO) and 30 meV (10Sc:BZO) for the 
quasielastic width. These values are considerably lower 
than the activation energy for the proton conductivity, 
which is around 500 meVP Thus, the process we ob- 
serve in our experiment cannot be rate limiting for the 
long-range proton diffusion. That is, there must exist 
another dynamical process of the protons, which occurs 
on a longer time scale and that cannot be observed in 
the present experiments, that controls the rate of the 
long-range diffusivity, i.e. the macroscopic proton con- 
ductivity. Perhaps, this could be the escape event from 
a trapping site in the structure, where the proton spends 
an extended time before it diffuses further. Furthermore, 
we note that the observed proton motion is slightly faster 
in 10Sc:BZO and that the quasielastic scattering inten- 
sity is slightly higher. Considering that 10Y:BZO has 
higher conductivity than 10Sc:BZO and a higher degree 
of hydration, 2 i.e. more protons in the structure, also 
point towards that we are observing a process that does 
not determine the rate of the long-range proton diffusion. 

As described above, we can imagine two localized dy- 
namical processes for the proton in hydrated perovskites, 
i) proton transfer between adjacent oxygens, and ii) ro- 
tation of the O-H group around the line connecting the 
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FIG. 6: Q-dependence of the quasielastic intensities at 380 K 
(o), 430 K (A), 465 K (0) and 495 K (□), respectively. Solid 
lines are fits to the data using a jump-diffusion model over 
two (transfer) and four (rotation) sites, according to Eq. |5| 
and Eq. (6}. 



octahedrally coordinated cations. Density functional the- 
ory (DFT) calculations have shown that these processes 
are well described with a jump diffusion model over two 
(transfer) and four (rotation) sites, respectively.^ For 
these models the quasielastic intensities are given by-ES 



1 sin(2Qr) 



and 



2 2 2Qr 



3 1 sin(2Qr) 1 sin(V2Qr) 



(5) 



2Qr 



V2Qr 



}, (6) 



where aL,a(Q) an d ol,4(Q) are the quasielastic structure 
factors for the transfer and rotational motion, respec- 
tively, r is the jump length, and a is a scale factor that 
we need to use since we do not have the quasielastic in- 
tensities on an absolute scale. Included in Fig. [6] are 
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free fits with the two model functions to the quasielas- 
tic intensity for the investigated temperatures, while the 
fit parameters are shown in Tab. |TJ As can be seen in 
Fig. [6] it is very difficult to discriminate between the two 
models. In order to directly separate them it would be 
necessary to have data at larger Q- values, ft should here 
be noted that some of us performed an experiment at the 
inverted time-of-flight quasielastic spectrometer IRIS at 
the pulsed neutron source ISIS in United Kingdom to ex- 
tend the Q _ran g e t° 0.5-3.7 A, but due to the presence 
of numerous Bragg peaks in this Q-range it was, unfortu- 
nately, impossible to analyze these spectra with sufficient 
accuracy. 

Insight into what process we are observing may instead 
be obtained from the jump lengths, see Tab. [I] From these 
we can see that for both the transfer and rotational model 
the jump lengths for the Sc-doped material are overall 
shorter than for the Y-doped counterpart. Furthermore, 
we see that apart from the lowest temperature (380 K), 
the obtained jump lengths are basically temperature in- 
dependent. The temperature independence for T = 430- 
495 K suggests that the fits are physically reasonable, 
while the overall larger values for the 380 K data is likely 
a consequence of the low quasielastic intensity at this 
temperature, wherefore these fit parameters have larger 
uncertainty. Thus, we will keep the following discussion 
to T > 430 K. 

The proton experiences a different local surrounding 
and hence a different potential energy surface (PES) de- 
pending on where it is located in the structure! 13 * 14 * Fig- 
ure [7] shows the relevant migration paths for the proton 
around and between oxygens in the first (01) and second 
(02) oxygen coordination shell of a dopant atom A. As 
seen in Fig. [7] Rl and R2 denote the rotational processes 
at 01 and 02, respectively, and Til, T12 and T22 rep- 
resent proton transfers between adjacent 01 atoms, be- 
tween adjacent 01 and 02 atoms, and between adjacent 
02 atoms, respectively. Further away from the dopant 
atom the PES is less distorted and there is essentially 
only one type of rotational (R) and one type of transfer 
(T) process. These local motions have been modeled by 
DFT calculations on a system corresponding to a dopant 
concentration of 3. 7% £21 The jump lengths and classical 
migration barriers so obtained for the Rl, R2, Til, T12, 
T22, R and T processes in Sc- and Y-doped BaZr03, re- 
spectively, are reported in Tab. [IT] By comparing these 
jump lengths to the experimental results in Tab. [TJ and 
by also taking the structure of the PES into account, 
we may draw conclusions about which process/es we are 
observing in the experiment. 

For the Sc-doped material the DFT calculations^ re- 
veal that the PES is much lower at and around 01, i.e. 
in the vicinity of the Sc atoms, than at 02 or at an O 
located even further away from the dopant. This means 
that the probability for finding a proton is higher close to 
a dopant than farther away, which in turn suggests that 
the local process observed in the experiment is either Rl 
or Til. The experimental jump length obtained from 



TABLE I: Jump lengths obtained from free fits of the 
quasielastic intensity with the jump-diffusion models over two 
and four sites according to Eq. ([5| and Eq. |6|. r2 denotes 
the jump length for the transfer process, and ta denotes the 
jump length (radius) for the rotational process. 



T/K 


10Sc:BZO 

r 2 /A r 4 /A 


10Y:BZO 

r 2 /A r 4 /A 


380 
430 
465 
495 


0.93±0.14 1.10±0.08 
0.83±0.10 0.98±0.06 
0.82±0.12 0.97±0.08 
0.82±0.10 0.96±0.08 


1.08±0.18 1.29±0.12 
0.99±0.12 1.18±0.10 
0.97±0.12 1.15±0.08 
0.97±0.12 1.14±0.08 



fitting with a transfer model, r 2 , is 0.82 A (c/. Tab. [i]), 
which is about 0.1 A longer than the DFT value for Til, 
given in Tab. nj The jump length (radius) obtained from 
fitting with the rotational model, — 0.97 A, is on the 
other hand in excellent agreement with the DFT value for 
Rl. From a structural point of view this would suggest 
that the observed local proton motion in 10Sc:BZO is the 
rotational motion around an oxygen close to a Sc atom. 
However, this conclusion is not supported by the classi- 
cal migration barriers (activation energies) presented in 
Tab.jlll The Rl process has a classical barrier of 210 meV, 
which is much higher than the 30 meV observed in the 
QENS measurements. Correcting for the difference in 
zero point energy, AZPE, of the proton at the stable site 
and at the barrier top is expected to reduce the calculated 
barrier with approximately 40 me Vf^ Consequently, this 
would still give a barrier of 170 meV, wherefore it is un- 
likely that Rl is the process we observe in the QENS 
experiments. 

For Til the DFT calculations give a much lower clas- 
sical energy barrier, only 130 meV compared to 210 meV 
for Rl. By further adding the zero point contribution 
(AZPE w -120 meV for proton transfers in BZ(j^) the 
theoretical barrier for Til is substantially reduced to ap- 
proximately 10 meV, which is close to the energy barrier 
determined from the experiments (30 meV). Based on 
this it is therefore more likely that it is the Til process 
that is observed. It should be noted, though, that the 
theoretical estimate most probably is an underestimation 
of the real Til barrier, as DFT at the level of complexity 
used in Re f.^ has a tendency to give too small barriers for 
hydrogen-bond mediated processes.^ Admittedly there 
is a discrepancy between the experimental and theoreti- 
cal jump lengths, see Tab. [T] and Tab. [Tl] but not so large 
that it cannot be explained by the errors inherent in the 
two methods. First of all, the experimental error bar is 
^0.1 A. Secondly, there is also an uncertainty in the cal- 
culated value. The ability of standard DFT approaches 
to accurately predict hydrogen bond lengths depends sen- 
sitively on the local geometry of the hydrogen bondP^ 
Often, but not always, the hydrogen bond length i s un- 
derestimated, sometimes with as much as ~ 0.1 

For the Y-doped material the situation seems to be 
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FIG. 7: Schematic picture of the migration paths for a proton 
in the vicinity of a dopant atom A, redrawn from RefP^l 

TABLE II: Jump lengths, r, and classical energy barriers 
(activation energies), E & , obtained from DFT calculations for 
the different migration paths in doped BZOP^Rl, . . ., R2 
are the paths in the near dopant region, as depicted in Fig. [7] 
whereas R and T are the rotation and transfer paths, respec- 
tively, far from dopant atoms. 



Migration 


3.7Sc:BZO 




3.7Y:BZO 




path 


r/A 


£ a /meV 


r/A 


£a/meV 


Rl 


0.96 


210 


0.97 


260 


Til 


0.72 


130 


0.80 


90 


T12 


0.89 


570 


0.79 


280 


T22 


0.89 


210 


1.04 


340 


R2 


0.89 


220 


0.89 


320 


R 


0.98 


200 


0.98 


200 


T 


0.82 


170 


0.82 


180 



very similar although for this material the PES is not 
only lowered at 01 but also, to an equal extent, at 02™ 
This suggests that the protons spend most of their time 
close to the Y dopants, but, in contrast to Sc:BZO, with 
a similar occupation of 01 and 02 sites. However, in- 
spection of Tab. [Tl] shows that all the migration path- 
ways Rl, . . ., R2 near Y, except Til, are associated with 
relatively high classical energy barriers (~ 300 meV). 
Hence, these processes are expected to occur on much 
longer time scales than the one/ones we observe in the 
experiments. The only remaining candidate is then Til, 
which has a calculated classical barrier is of 90 meV (c/. 
Tab. 0. Further addition of AZPE ps -120 meV obvi- 
ously results in a negative barrier, which is not physically 
sound. Given the tendency of DFT to underestimate this 
type of barrier it is, however, not unreasonable to link 
Til to the measured 10 meV activation energy in the 
Y-doped material. Consequently, we once again claim 
that Til gives the predominant contribution to the ex- 
perimental spectra, despite an even larger (0.17 A) dif- 
ference between the experimental and theoretical jump 



lengths. In this context we note that even though the 
degree of hydration is higher in the Y-doped materialpl 
the slightly higher quasielastic intensity for the Sc-doped 
material (c.f. Fig. [6]) suggests that there is a larger num- 
ber of protons in the immediate vicinity of the dopant in 
the Sc-doped material than in the Y-doped counterpart. 
This is in agreement with the very high T12 barrier for 
the Sc-doped material, which may also be the reason to 
the lower proton conductivity of this material. 

Finally, it should be noted that the comparatively 
larger disagreement between our experimental results 
and the calculations in 10Y:BZO may suggest that a 
more complex process is taking place, or that the pro- 
cesses listed in Tab. |TT| contribute in an averaged manner 
to the quasielastic intensity. In order to elucidate what 
process/es we are seeing in the two materials, and in 
10Y:BZO in particular, further experimental and theo- 
retical investigations are needed. 

VI. CONCLUSIONS 

Proton motions in the hydrated proton conducting per- 
ovskites BaZro.9Ao.10O2.95 (A — Y and Sc) have been 
investigated using quasielastic neutron scattering. The 
results reveal a localized proton motion on the ps time 
scale and with a low activation energy, ~10-30 meV. 
Comparison of the QENS results to density functional 
theory calculations suggests that for both materials we 
observe proton transfers between neighboring oxygens 
bound to a dopant atom. The low activation energy, 
much lower than the activation energy for the macro- 
scopic proton conductivity, suggests that this process is 
not rate-limiting for the long-range proton diffusion in 
either of the two materials. That is, there must exist an- 
other dynamical process of the protons, which occurs on 
a slower time scale and that cannot be observed in the 
present experiment, that controls the rate of the long- 
range diffusivity, perhaps the escape event from a trap, 
where the proton spends a longer time before it diffuses 
further. So, even though the two investigated materi- 
als possess considerably different proton conductivities, 
this difference cannot be directly linked to the observed 
proton motion. 
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